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mOBAJBLE SOUR FIARE DOSES ENCOUNTERED ON AN INTERPIANETARY MISSION AS GALCUUTED BY THE MCFURE CODE 


Gerald P. Lahti and Irving M. Karp 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 

The computer program, MCFURE, uses Monte Carlo methods to simulate solar 
flare occurrences during an interplanetary space voyage. The total biologic al 
dose inside a shielded crew compartment due to the flares encountered during 
the voyage is determined. The computer program evaluates the doses obtained 
on a large n\amber of trips having identical trajectories. From these results, 
a dose Dp having a probability p of not being exceeded during the voyage 
can be determined as a function of p for any shield material configuration. 

The user of the code selects any number of solar flares considered to be 
representative of the ones that will occur during future solar active periods 
(these flares are generally selected from the flares that occurred during the 
last solar active period 1956 to 1962). The flares are assumed to occur dur- 
ing these periods. The dose at a distance of 1 AU from the Sun from each of 
these flares behind any shield configuration investigated is input to the 
MCFIARE code. The code accounts for the dependence of the dose received from 
a flare on the distance from the Sun according to a (l/r)°' variation, where 
r is the distance from the Sun and the exponent a can be assigned. From 
trajectory parameters, which are input to the computer program, the distance 
from the Sun as a function of time during the trip is calculated. 

To illustrate the use of the code, a trip to Mars and return is calcu- 
lated, and estimated doses behind several thicknesses of aluminum shield and 
water shield are presented. 


INTRODUCTION 

The protons emitted by solar flares are the 
most hazardous source of space radiation encount- 
ered on an interplanetary space voyage which rapid- 
ly traverses the Van Allen belts. To define shield- 
ing requirements for biological protection from 
solar flare protons encountered in manned inter- 
planetary missions, a computer program MCFIARE 
(ref. l) has been prepared which uses Monte Carlo 
methods to simulate solar flare occurrences during 
the trip and records the toteil associated biological 
doses for the trip behind various shield material 
configurations. 

A set of flare events, considered to be repre- 
sentative of those that will occur during future 
solar active periods, is required. Presently, this 
set is selected from the flares that occurred dur- 
ing the last solar active period from 1956 to 1962. 
The flares are assumed to occur randomly during the 
trip. The doses behind the various shield config- 
urations investigated have to be evaluated for each 
flare, at a distance of 1 AU from the Sun, by using 
another computer code such as one of those described 
in Refs. 2 or 3. (Ref. 4 presents doses from solar 
flare proton spectra behind two shield materials as 
calculated by using Ref. 2.) These dose values are 
input to the computer program MCFIARE. This pro- 
gram provides for including the dependence of the 
magnitude of the dose received from the flare on 
the distance r from the Sun. This dependence is 
taken to be of the form (rQ/r)°^, where Tq is the 
distance between the Sun and the Earth, and the ex- 
ponent a can be any arbitrarily selected value. 

The Monte Carlo sampling procedure (analogous 
to particle transport Monte Carlo) is used to de- 
termine when flares are encountered on the trip, 
which of the input flares occur, and the total dose 
obtained from all flares encountered during the 
trip behind each shield configuration. The comput- 
er program tallies the doses obtained for a large 
member of trips having identical trajectories and 
from this, the distribution of doses incurred on 
the mission for each shield configuration is ob- 


tained. Each dose distribution determines a dose 
Dp having the probability p of not being exceed- 
ed as a function of p. 

METHOD OF ANALYSIS 

Any N flare events can be selected as repre- 
sentative of the flare activity that will occur 
during the space voyage. The doses from each flare 
at a distance of 1 AU from the Sun behind one or 
more shield material configurations investigated 
are input to the computer. Each flare is assumed 
to occur randomly on the average of once during the 
solar active period T^gf, which is of the order of 
2000 days or 5^ years. 

Monte Carlo Method to Determine Occurrences 
of Flares During Trip 

Inasmuch as flares are assumed to occur ran- 
domly in the solar active period, the probability 
of any flare occurring in the time interval dt is 
(N/r^ef H dt. Then the probability of no 
flare occurring for time t is e“^A■t and the 
probability of a flare occurring in dt at time t 
is dt. 

In the Monte Carlo method, this probability of 
flare occurrence in the time interval dt at time 
t is equated to the probability of selecting a 
random number in an interval d^ about If a 

random number set uniformly distributed in the in- 
terval (0,l) is selected, then the probability of 
selecting a random number in d^ about ^ is d^. 
Equating these two probabilities, one gets 

dt = d? 



t = -1 ln(l - 0 


1 



( 3 ) 


Inasmuch as both | and 1 - | are uniformly 
distributed in the interval (0,l) an equivalent ex- 
pression for t which results in the same exponen- 
tial distribution is 




/p(i - p) 

V Mf2(Dp) 


t = - 1 In I (1) 

By selecting a random niomber from a uniformly dis- 
tributed set in the interval (0,l) the time be- 
tween flare events is obtained from Eq. (l). 

When an event has occurred, a’ particular flare 
of one of the N is also selected by random num- 
bers. The i"^^ flare occurs when a new random 
number occurs in the interval 

x = l, 2, . . N (2) 

Dose values input to the code are calculated 
by using the methods of Ref. 2, or 3, for each rep- 
resentative flare using the solar flare proton in- 
tensity and spectral shape as observed at the 
Earth. The code permits an effect of the distance 
from the Sun on the dose (or intensity) received 
from a flare of the form where r^ is the 

distance from the Sun to the Earth, r is the dis- 
tance of the spacecraft from the Sun during the 
flare, and a is an exponent that is assigned. 

The code evaluates r for any time during the 
trip by calculating the mission trajectory from in- 
put parameters describing the two transfer ellipses 
and the total trip time. The first ellipse is the 
trajectory from the Earth to the planet and the 
second describes the return trip. 

The calculation proceeds as follows: Select- 

ing a random number and using Eq. (l) determines 
the time in the mission when a flare occurs. Se- 
lecting another random number and using Eq. (2) de- 
termines the particular flare which occurs. Know- 
ing the time in the mission (and, hence, distance 
from the Sun) and the particular flare, the dose 
received through each shield configuration is 
tallied. The duration of the flare event is then 
added to the time in the mission when the flare 
occurred to give the time when the event is over. 
From a new random mmiber and Eq. (l), the time 
elapsed until the next flare occurs is calculated, 
and so on. This procedure is repeated until the 
total trip time has elapsed. The total dose from 
all flare events encountered on the trip is deter- 
mined for each of the shields considered. 

Selection of Dp 

This procedure is repeated for a large number 
M of trips of identical trajectory. The code di- 
vides the dose range into dose intervals that are 
1 dose unit wide and tabulates the number of trips 
A(D) during which doses that lie within each dose 
interval are encountered. The cumulative fraction 
of the total trips that encounter doses less than 
the upper dose bound of each interval is also tab- 
ulated (this is an estimate of the probability p 
that a dose equal to the upper dose bound will not 
be exceeded). From this tabulation, for any spec- 
ified p, a dose Dp which has the probability p 
of not being exceeded on the trip can be selected 
for any shield. The standard deviation of this 
value from the true value is shown in Ref. 1 to be 


where M is the number of trips in the group in- 
vestigated, and f(Dp) is the dose probability 
density at Dp (which is calculated from the tab- 
ulated output). If A(Dp) is the number of trips 
contained in the interval in which Dp occurs, 
f(Dp) is approximately given by 


f(Dp) 


A(Dp) 

M 


If A(D) is a widely fluctuating value in the vi- 
cinity of Dp, averaged values of A(D) can be hand 
calculated from the tabular computer output; 
f(Dp) is obtained from the averaged value of 
A(D ). Note that for a given p, because f(Dp) 
is independent of the number of trips, varies 

inversely as y/u. 

To reduce the uncertainty associated with Dp 
for a given number of trips, the code is written 
so that the trips constitute a stratified sample 
for the time-to-first-flare selected in accordance 
with the exponential distribution exp[-(Nt/Tref)J • 

ILLUSTRATIVE EXAMPLE 

An interplanetary trip to Mars and return has 
been evaluated in order to illustrate details of 
the method and to present some representative 
shield requirements for such a mission. Water and 
aluminum were the shield materials considered. 

Selection of Flare Events 


The flare occurrences during the last solar 
active period 1956 to 1962 are assumed to be repre- 
sentative of those that will occur during a future 
active period. Webber in Ref. 5 has compiled a 
record of these flares, their time -integrated 
proton intensities, and spectra. From these rec- 
ords, the 20 largest flares, based on proton inten- 
sity, were selected (the effect of neglecting the 
rest of the flares is small for shields considered 
here). These flares are listed in table I which 
presents the integrated flux of protons having 
energies greater than 30 MeV and some constants A 
and pQ associated with each flare. The spectral 
shape of each flare is assimied to vary as 


N > E == A exp 


(- 




+ 1876E 


where N > E is the total number of protons per 
square centimeter with energies greater than E. 

The values of Pq were obtained from Webber's 
compilation, and the values of A were calculated 
to be consistent with flux values of N > 30 MeV. 

It has been observed that some of the flares 
tend to occur in clusters. Table I indicates that 
such clusters occurred in August 1958, July 1959, 
November 1960, and July 1961, These four clusters 
were selected as representative of the clustered 
flare events that may occur and were assumed to be 
events of 12-, 8-, 10 -, and 8 -day durations, re- 
spectively. The other nine single flares were each 
assumed to have a duration of 2 days. During fu- 
ture solar active periods, each of the four clust- 
ered events and each of the nine single flare 
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events are assumed to occur randomly with a fre- 
quency of once during the active period (taken to 
he 2000 days). 

Doses For Selected Flare Events 

Table II presents the doses obtained behind 
various shield thicknesses from each of these 13 
events. These doses were calculated by using the 
Lewis Proton Shielding Code^ described in Ref. 2 , 
and are representative of doses received at the 
center of a spherical crew compartment having the 
given shield thicknesses when the vehicle is lo- 
cated at a distance of 1 AU from the Sun. 

Table Il(a) presents the rem doses behind the water 
shield^ and table Il(b) the rem doses behind the 
aluminum shield. The data presented in this table 
are input to the computer code MCFLARE. 

Trajectory 

The trajectory selected for the trip is one 
that requires 556 days for the complete mission 
(and happens to be one that results in near mini- 
mum vehicle weight for a departure date in 1983) . 
The outward journey from Earth to Mars requires 
280 days, then there is a 40-day stay at Mars, then 
a 236-day return trip during which the vehicle ap- 
proaches within 0.5 AU of the Sun. 

Discussion of Results 

Figure 1 shows the distribution of dose ob- 
tained for a large number of identical trips. The 
figure presents the rem dose distribution behind 
20 g per square cm of water shield obtained from a 
computer run of 40 000 trips. The effect of dis- 
tance on dose was assumed to vary as (l/r)^. 

In Fig. 1(a), A(D), the number of trips on 
which a dose between D - 1 and D has been en- 
countered, is plotted against D. A smoothed 
curve of this distribution A(D) (obtained by av- 
eraging seven values centered about each A(D)) is 
also shovna as a dashed line in this figure. 

In Fig. l(b), p(D) the probability of not ex- 
ceeding any dose D is shown plotted against D 

p(D) = 1 

i=l 


where M is the number of trips considered in a 
computer run. For the case of Fig. l(b). 


P(D) 


40 000 


D 

£ A(i) 


i=l 


If Dp is defined as that dose which has a 
probability p of not being exceeded, then from 
Fig. l(b), Dq .99 for the group of 40 000 trips, 
occurs in the dose interval between 93 and 94 rem 
and is selected as the upper bound of the interval, 
namely, 94 rem. 

The dose probability density function f(D) is 
approximately equal to A(D) /M. At a Dq. 99 of 
94 rem, the value of f(D), evaluated by using the 
smoothed value of A(D), is 


and the standard deviation from Eq. (3) is 

^Do.99 " 

There is a 97.7 percent confidence that the 
true dose Dq, 99 is less than D 0.99 + 2 opQ^gg; 
is less than 95.5 rem. 

Table III presents values of D 0.99 deter- 
mined by using the computer code for rem doses be- 
hind various thicknesses (expressed as g/cm^) of 
aluminum shield and water shield. These values 
are presented for both the case where the dose re- 
ceived from a flare event is assimied to vary as 
(l/r)^ (considered representative of what might be 
a solar proton diffusion model in space) and for 
the case where there is no effect of distance 
(i.e., distance effect, g(r) = l) . Also shown in 
the table are the values of associated 

with each Dq .99 corresponding values of 

1^0.99 + 99- 

In Fig. 2 the values of Dq. 99 + 2 aDo, 9 Q 
plotted against shield thickness. In Fig. 2(a), 
the rem dose is plotted against water shield thick- 
ness for cases of g(r) = l/r^ and g(r) = 1 . 
Figure 2(b) is a similar plot of the rem dose be- 
hind the aluminum shield. From curves such as 
Pig. 2, one can select a shield thickness such that 
on a given mission, there is a probability p of 
not exceeding an accumulated dose Dp. 

These curves indicate how much more effective 
water is as a shield material than aluminum for 
shielding against solar flare protons on the basis 
of grams per square centimeter of material neces- 
sary to maintain a given dose level inside the 
crew compartment. Also shown for this mission is 
that a l/r^ distance effect on dose received 
from a flare can have an appreciable effect on 
shield requirements when the mission trajectory 
brings the vehicle to within 0.5 AU of the Sun. 

This result indicates a need for information re- 
garding the effect of position from the Sun on ra- 
diation encountered from a flare, 

COMPUTER PROGRAM 

Complete data input instructions for the com- 
puter program MCFLARE are presented in Ref. 1. 

The version of the code presented is operational 
on the Lewis Research Center IBM 7094-11/7044 com- 
puter system. Execution times for the examples 
discussed previously in this report were about 
2 minutes per group of 40 000 trips. The MCFLARE 
code is available from the Radiation Shielding In- 
formation Center of Oak Ridge National Laboratory 
as code package CCC-93. 


f(Do 99 ) = 26 = 0.00065/rem 

40 000 
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TABLE I. - FLARES (FROM SOLAR ACTIVE PERIOD 1956 TO 1962) 
SELECTED AS REPRESENTATIVE OF THOSE THAT WILL 


OCCUR DURING FUTURE ACTIVE PERIODS 

[Integral flux above 30 MeV, and spectral constants A and P^ for 
each flare are tabulated. ] 


Flare 

date 

N > 30 MeV 
protons/ cm^ 

^0 

A, 

protons/ cm^ 

2-23-56 

l.OKlO^ 

195 

3.4x10® 


1-20-57 


61 



8-29-57 

1.2x10^ 

56 

8, 6x10^ 



5x10^ 

127 

3. 3x10® 


3-23-58 

2,5x10® 

64 

1. 1x10® 


7-7-58 

2, 5x10® 

62 

1.2x10® 


8-16-58 

4x10^^ 

64 

1.7x10® ^ 


8-22-58 

Txio"^ 

56 

5,0x10^ 

► Clustered flare event 

8-26-58 

1. IxlO® 

51 

1.2x10^^ 



9.6X10® 

84 



6-13-59 

8. Sxio"^ 

^48 



7-10-59 

1,0x10® 



1 

7-14-59 

1.3x10® 

80 


^ Clustered flare event 

7-16-59 

9, lx 10® 


8.9x10^ ^ 

1 

9-3-60 

3. 5Xlo'^ 

127 

2.3x10® 


11-12-60 

1.3X10® 

124 

8. 9x10® 


11-15-60 

7.2x10® 

114 

q 

5.9X10^ 

► Clustered flare event 

11-20-60 

4, 5xlo’ 

118 

3.4x10® ^ 


7-12-61 

4xlo'^ 

56 

2,9x10® 

^ Clustered flare event 

7-18-61 

3x10^ 

102 

3. 1x10^ J 

[ 


^Estimated, 






TABLE n. - REM DOSES FROM SELECTED FLARE EVENTS BEHIND 


VARIOUS THICKNESSES OF WATER AND ALUMINUM SHIELD AT 
1 ASTRONOMICAL UNIT FROM SUN 


(a) Rem dose behind water shield 


Shield thickness, 
g/cm^ 

Flare date 

10 

15 

20 

30 

40 

Dose, 

rem 


27.76 

20. 79 

13.03 

9. 03 

2-23-56 


.38 

.19 

.07 

.04 

1-20-57 


.15 

.07 

.03 

.01 

8-29-57 


.65 

.42 

.21 

.12 

10-20-57 


.40 

.19 

.07 

.04 

3-23-58 

1.26 

.45 

.22 

.08 

.04 

7-7-58 

1.03 

.35 

.17 

.07 

.04 

Aug. 1958 cluster 

11.30 

4.93 

2.67 i 

1.09 

.57 1 

5-10-59 

.60 

.20 

.10 

.04 

.02 

6-13-59 

49.59 

23.68 

13.76 

6.11 

3.36 

July 1959 cluster 

.84 

.46 

.29 

.15 

.09 

9-3-60 

49.34 

26.49 

16.79 

8.30 

4.85 

Nov. 1960 cluster 

5.48 

2.67 

1. 57 

.71 

,39 

July 1961 cluster 


10 


53.93 
2.11 
.90 
1.73 
2.20 
2.48 
2.12 
18.98 
1.25 
78.42 
1.22 
71.89 
8. 57 


(b) Rem dose behind aluminum shield 


Shield thickness, 

g/cm^ j 

Flare date 

15 



20 

30 

40 

50 

60 


Dose, 

rem 

37.20 

28.77 

19.07 

14. 16 

11.15 

9.26 

2.23-56 

.86 

.48 

.24 

.16 

.12 

.10 

1-20-57 

.35 

.19 

.10 

.07 

.06 

.05 

8-29-57 

1.00 

.68 

.38 

.25 

.18 

.14 

10-20-57 

.89 

.50 

.25 

.17 

.12 

.11 

3-23-58 

1.01 

.56 

.28 

.19 

.14 

.12 

7-7-58 

.83 

.45 

.24 

.17 

.13 

.11 

Aug. 1958 cluster 

9.06 

5.36 

2.59 

1.63 

1.18 

.93 

5-10-59 

.49 

.26 

.14 

.10 

.08 

.06 

6-13-59 

40. 39 

25. 19 

12.79 

8. 12 

5.86 

4. 59 

July 1959 cluster 

.71 

.48 

.26 

.17 

.13 

. 10 

9-3-60 

41.17 

27,59 

15.09 

9.87 

7.17 

5.61 

Nov. 1960 cluster 

4.49 

2.82 

1.44 

.92 

.66 

.52 

July 1961 cluster 






























TABLE m. - VALUES OF Dq gg AND CORRESPONDING STANDARD 
DEVIATION BEHIND VARIOUS THICKNESSES OF WATER AND 
ALUMINUM DETERMINED FROM GROUP OF 40 000 TRIPS 
(a) Rem dose behind water shield 


Distance 

effect, 

Quantity 


Shield thickness, 
g/cir? 


g(r) 


10 

15 

20 

30 

40 





Dose, 

rem 



l/r2 

°0.99 

254 

139 

1 

94 

55 

37 


X99 

2.0 

1.3 

.77 

.44 

.31 


*^0.99 + 2^^0 99 

258 

141.6 

95. 5 

55.9 

37.6 

1 

°0. 99 

165 

90 

63 

35 

23 


X99 

1.7 

1.0 

.57 

.25 

.2 


°0.99 + 2<TDooo 

168.4 

92 

64. 1 

35. 5 

23.4 


(b) Rem dose behind aluminum shield 


Distance 

Quantity 



Shield thickness. 



effect, 





g/cm^ 




g(r) 


10 

15 

20 

30 

40 

50 

60 






Dose, 









rem 





%99 

376 

213 

146 

87 

62 

48 

39 


99 

3.3 

1.8 

1.4 

.77 

. 54 

.42 

.31 


Do. 99 +2^00 09 

382.6 

216.3 

148.8 

88. 5 

63. 1 

48.8 

39.6 

1 

°0.99 

250 

m 

94 

58 

40 

31 

25 


X99 

2.2 

H 

.91 

.59 

.31 

.23 

.22 


D0.99+2<^Dg gg 

254.4 

m 

95.8 

59.2 

40.6 

31.5 

25.4 
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CUMUUmVE PROBABILITY, p(D) 


1.00 
.98 
.96 
.94 
.92 
.90 
.88 
.86 

30 40 50 60 70 80 90 100 110 120 130 

DOSE, D, REM 

(B) CUMULATIVE PROBABILITY OF TRIP ENCOUNTERING DOSE LESS THAN D. 

Figure 1. - Concluded. 




Figure 2. - Values of Dq 99 + 2aQ^ for various 

thicknesses of water and aluminum shields, and 
for cases of distance effect varying as 1/r^ and 
no effect of distance. Number of trips, 40 000. 



DOSE, D, REM 



(B) REM DOSE BEHIND ALUMINUM SHIELD. 
Figure 2. -Concluded. 


